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In  this  paper  we  continue  the  study  of  balanced  treatment  incomplete 


block  (3TI8)  designs  which  wu  initiated  in  [3}\  These  designs  are 
appropriate  for  comparing  simultaneously  p  ">  2  test  treatments  with 

~y  r>> 

a  control  treatment--the  so-called  multiple  comparisons  with  a  control 
(MCC)  problem.  This  class  of  designs  was  characterized  in  [3],  In  the 
present  paper  we  obtain  optifnal  designs  within  this  class  for  selected 
(p,k,b)  where  k  <  p+1  is  the  number  of  plots  per  block,  and  b  is 
the  total  number  of  blocks  available.  Specifically,  optimal  designs 
are  obtained  for  d  =  2(1)6,  k=2  and  for  p  =  3,  k  =  3. 

Tables  of  exact  (discrete)  optimal  designs  are  given  for  these 
(p,k)-values  for  a  range  of  b-values  which  would  ordinarily  be  of 
practical  interest.  Tables  of  approximate  (continuous)  optimal  designs 
are  given  for  situations  in  which  very  large  b-values  are  required.  The 
theory  underlying  these  approximate  designs  is  developed,  and  the  goodness 
of  the  approximation  is  studied. 

- - 

Key  words  and  phrases :  Multiple  comparisons  with  a  control,  balanced 
treatment  incomplete  block  (BTIB)  designs,  admissible  designs,  optimal 
designs,  applications  of  equicorrelated  multivariate  normal  and  multi¬ 


variate  Student's  t  distributions. 


1.  INTRODUCTION  AND  SUMMARY 


In  our  first  paper  [3]  we  initiated  the  study  of  balanced  treatment 
incomplete  block  (BTIB)  designs  which  are  appropriate  for  comparing 
simultaneously  p  >  ^  test  treatments  with  a  control  treatment.  This 
class  of  designs  was  characterized  in  [3];  in  the  present  paper  we 
obtain  optimal  designs  within  this  class  for  selected  (p,k,b)  where 
k  <  p+1  is  the  number  of  plots  per  block,  and  b  is  the  total  number 
of  blocks  available  for  experimentation. 

In  order  to  make  this  paper  self-contained  we  will  state  the  key 
definitions  and  results  from  [3],  We  shall  use  the  following  notation 
(also  used  in  [3]):  Let  the  treatments  be  indexed  by  0,l,...,p  with 
0  denoting  the  control  treatment  and  1,0,..., p  denoting  the  test 
treatments.  The  N  =  kb  experimental  units  can  be  arranged  in  b 
blocks  each  of  size  k.  If  treatment  i  is  assigned  to  the  hth  plot 
of  the  jth  block  (O^i^p,  H  h  <  kj  1  <  j  <  b),  let  Y .  denote 

the  corresponding  random  variable;  we  assume  the  usual  additive  linear 
model  (no  treatment'^block  interaction) 


Yi)h =  “  *  *  ej  *  nih  u-n 

p  b  ^ 

with  ^  a.  =  7  8.  =  0;  the  e...  are  assumed  to  be  i.i.d.  N(0,o‘) 

i*0  1  j-1  j  .  ljh 

random  variables.  It  is  desired  to  make  an  exact  joint  confidence  state 
ment  (employing  one-sided  or  two-sided  intervals)  concerning  the  p 
differences  a,  -  a.  based  on  their  best  linear  unbiased  estimators 

0  i 

(BLUE’s)  a.  -a,  (1  <  i  <  p). 


In  Section  3.1  of  [3]  we  proposed  a  class  of  incomplete  block  designs 

which  are  balanced  with  respect  to  the  test  treatments  in  the  following 

-  .  2  2  * 
sense:  Varla,  -  a.)  s  n  o  /N  (1  '  i  <  p)  and  Corria,  -  a.  ,  a.  -  a.  I 
Oi  ~  “*  0  i.  0  i- 

(i,  *  i,;  1  <  i,  ,i,  <,  p);  the  parameters  n  and  o  depend  on  the  design 
employed.  We  refer  to  designs  with  this  property  as  3TI3  designs.  Con¬ 


ditions  that  a  design  must  satisfy  in  order  that  it  be  BTI3  were  given 


in  Theorem  3.1  of  [3].  This  theorem  states  that  if  {r.^}  is  the 

incidence  matrix  of  the  design,  r. .  being  the  total  number  of  times 

11  b 

the  ith  treatment  appears  in  the  jth  block,  and  if  X.  .  =  £  r.  .r.  . 

1112  «  =  1  Al3  i2" 


is  the  total  number  of  times  that  the  i^th  treatment  appears  with  the 
i.th  treatment  in  the  same  block  over  the  whole  design 


ii,  *  i,;  0  i  ,i,  ^  p),  then  the  necessary’  and  sufficient  conditions 

L  «.  X  A 

for  a  design  to  be  BTI3  are  that 


In  Section  »*  of  [3]  we  restricted  consideration  to  3TI3  designs,  and  shewed 
how  to  use  such  designs  for  experiments  leading  to  joint  one-sided  (.or 
two-sided)  confidence  interval  estimates  of  the  a  -  a.  (1  ^  i  <  p) 

U  1. 

when  o“  is  known  or  unknown. 

The  concept  of  an  admissible  design  (see  Section  3.1  below)  for  such 
experiments  was  introduced  in  Section  5  of  [3].  The  problem  of  finding 
an  optimal  design  in  the  class  of  admissible  BTI3  designs  when  ip,k,b) 
are  given  was  mentioned  in  Section  <*  of  [3]  but  was  not  precisely  formu¬ 
lated.  We  do  so  in  Section  3.2  of  the  present  paper,  and  solve  the  problem 


3 


for  the  cases  p  =  3(1)6,  k  =  3  and  p  =  3, 
combinations  of  practical  interest  with  k  ^  3 
The  specific  MCC  problem  with  which  we  are 
paper  is  that  of  obtaining  point  one-sided  cent 


k  =  3;  additional  tp,k>- 
wili  be  considered  in  [4]. 
concerned  in  the  present 


idence  intervals  or  the 


form 


a  -  a.  >  a, 
3  1  —  0 


a .  -  d  ( 1  <  i  <  p ) 
1  -■  —  ' 


(1.3) 


for  given  values  of  (p,k,b)  when  0*  is  known,  and 


is  a  speci¬ 


fied  "vardstick"  associated  with  the  "width"  of  the  confidence  interval. 


For  this  problem  we  seek  an  opt imal  design  in  the  class  of  all  admissible 
5715  designs,  an  optimal  design  being  one  which  maxim  ices  the  confidence 
coefficient  associated  with  tl.3). 

3.  AN  INTEGRAL  EXPRESSION  FOR  THU  C0NFIDCN3E  COEFFICIENT 

For  ease  of  reference  we  record  here  the  expressions  derived  in  [3] 

for  the  estimators  a  -  a.  (1  <  i  <  p),  and  their  variances  and  cor- 

relations.  Let  T.  denote  the  sum  of  all  observations  obtained  with 
1 

the  ith  treatment  (0  <  i  <  p),  and  let  5.  denote  the  sum  of  all 

5  b 

observations  in  the  ith  block  (1  *■  i  <  b).  Define  5*?  =  ^  r..5. 

”  ■  “  1  .  1 '  - 

*  3  =  - 

and  let  Q.  =  kT.  -  3.  (0  <  i  <  p).  Then 

1  1  1  —  -  ' 


xiQo  -  V\  ...  , 

ao  ■  ai  -  \,r\  tr\,)  a  i 1  - 
0  01 


r  3 . 1 ) 


Also , 


Var{a,  -  a.)  =  ^r-  o'  (1  v  i  <  p) 
U  1  N  *  • 


(3.3) 


4 


where 


ar. 


n 


V  Vrvi’  • 


0 


Corr(a. 

u 


l  i,  *  i 


p).  (2.4) 


The  probability  associated  with  (1.2)  is  given  by 


where 


P(a ,  -  a ,  >  i 


j  i  ai3  ■  ^  ( 1  £  i  <,  p)} 


=  P{Zi  <  d.'N  on  (1  <  i  <  p)) 


12.5) 


(a, -a .  )  -  (a  -a .  ) 


0  i. 


on/*  N 


(1  ^  i  i  p) . 


The  2^  have  a  standard  p-variate  normal  distribution  with 

Corr{^.  ,w.  }  =  o  (i^  *  ij)  i  i  i^.i-,  2.  ?);  n*  and  p  are  given  bv 
L  *. 

(2.3)  and  (2.4),  respectively.  We  next  define 


5  =  d*M/o 


(2.6) 


which  is  a  pure  number  involving  the  specified  quantity  d.  Then  (2.5) 


can  be  written  as 


'lx) 


(2. 


) 


where  $( • )  represents  the  standard  normal  cdf. 


FORMULATION  OF  THE  EXACT  (DISCRETE)  OPTIMIZATION  PROBLEM 


3.1  Inadmissible,  admissible,  equivalent  and  generator  designs 

Preliminary  ro  posing  our  optimization  problem  we  recall  some 

definitions  from  [3].  If  for  given  (p,k,b)  we  have  two  3TI3 

2  2 

designs  D,  and  D.  with  (n,,P,)  and  (n-.p-),  rescectively , 

X  7  1  i  1 

2  2 

and  if  <,  and  ^  c0  with  at  least  one  inequality  being 
strict,  then  D„  is  said  to  be  inadmissible  since  it  is 
dominated  (in  the  sense  of  having  a  smaller  confidence  coefficient 
for  every  d  and  a)  by  D  ;  if  a  design  is  not  inadmissible, 

then  it  is  said  to  be  admissible.  (A  more  general  definition  of  inadmis- 

2  ° 

sibility  is  given  in  [**].  See  also  Remark  3.1,  below.)  If  n,  = 

1 

and  p,  =  o„  then  D,  and  D„  are  said  to  be  equivalent;  if 

-  l  7  — * — 

(\!)")  ,11*))  are  the  parameters  associated  with  (i  =  1,2),  then 

D,  and  D,  are  equivalent  if  and  only  if  and  =  \^x^. 

To  construct  3TI3  designs  for  given  (p,k)  and  any  b  we  next 
introduce  the  concept  of  a  generator  design.  For  given  (p,k)  a 
generator  design  is  a  3TI3  design  no  proper  subset  of  whose  blocks 
forms  a  3TI3  design.  The  candidates  for  an  optimal  design  for  given 


(p,k,b)  and  specified  d/o  will  be  ail  of  the  admissible  3T3 


2  S  1 


that  can  be  constructed  for  given  b  from  the  list  of  all  generator 
designs  for  given  (p,k).  If  two  or  mere  equivalent  admissible  3TI3 
are  available,  then  it  is  necessary  to  consider  only  one  of  them  (for 
the  purposes  of  solving  the  optimization  problem  (3.1),  below);  thus 


in  general  for  given  (p,k)  we  need  to  consider  only  the  nonequiva- 
lent  generator  designs  in  our  search  for  an  optimal  design.  For 
example,  for  (p,k)  =  (u,3)  the  designs 


6  a  /b 


with  X^1*  =  X(‘J  =  2,  XU)  =  =  2  are  equivalent 

0  0  1  - 

designs,  but  only  one  of  them  need  be  considered  in  the 
problem.  As  a  second  example,  for  (p,k)  =  (4,3)  the 


generator 

optimization 

design 


is  a  generator  design  with  X  ^  =2,  X,  '  =  2  which  is  equivalent  to 

the  union  of  the  two  generator  designs 


with  xi4)  =  2,  x|4)  =  0,  X<5)  =  0,  X^5)  =  2,  respectively,  thus 

0  10  1 

it  suffices  to  consider  the  generator  designs  D,,  and  D  in  the 

4  D 

optimization  problem. 


MMMNMMmm 


6  b/b 


/  1  \  (  '  ) 

Remark  3.1:  We  note  that  (X„  , X,  >  =  (2,2)  and 

-  -  ~  0  1 

(\!;4)  +  X^},  x!'*)  +  X^:))  =  (2,2)  for  the  3TIB  designs  2,  and 


respectively.  Hence  7ar{ 


‘.cu  -  a.  ; 


;f  (2.2)  is  the  same 


for  2,  and  3  u  D_  (since  N  =  kh  and  hence  n  /N  in  (2.2)  is 


*  4  3 

the  same  for  both  designs);  also,  p  of  (2.4)  is  the  sane  for  both 
designs.  This  implies  that  the  probability  (2.7)  is  the  same  for 
both  designs  even  though  D,  reauires  b  =  7  and  2  u  D„  recuires 

i  4  D 

b  =  8.  We  can  therefore  regard  D,  as  dominating  D4  u  Dc  since 
2,  recuires  a  smaller  number  of  blocks  than  2  u  D.  to  achieve 

1  •  -  4  3 

the  same  orobabilitv ;  in  this  sense  2,  u  D,  is  inadmissible  with 

-  -  4  3  -  - - 

respect  to  D  .  This  more  general  concept  of  inadmissibility  which 
allows  the  comparison  of  designs  with  different  b-values  is  introduced 
and  discussed  in  detail  in  C-O.  The  possibility  of  a  design  with  a 
smaller  b-value  dominating  a  design  with  a  larger  b-value  does  not 
arise  for  the  cases  p  =  2(1)6,  k  =  2  and  p  =  3,  k  =  3  considered 
in  detail  in  the  present  paper.  The  concept  of  strong  inadmissibility 
introduced  in  the  following  paragraph  also  is  generalized  in  [4]  to 
cover  situations  in  which  3TIB  designs  with  different  b-values  are  being 
compared. 


There  are  certain  generator  designs  which  will  always  yield  inadmis¬ 
sible  BTI3  designs,  and  we  would  like  to  eliminate  such  designs  in  our 
search  for  an  optimal  design;  we  call  such  designs  strongly  (S-)  inadmis¬ 
sible  designs.  If  for  given  (?,k)  we  have  two  designs  ,  D2  (not 


7 


necessarily  generator  designs)  with  b,  =  b2 ,  we  say  that  D0  is 
S- inadmissible  with  respect  to  D  if  0o  is  inadmissible  with  respect 
to  D^,  and  for  any  arbitrary  BTIB  design  we  have  that  u 

is  inadmissible  with  respect  to  u  D^.  S-inadmissibility  implies 
inadmissibility  but  not  conversely,  for  example,  for  (p,k)  =  (4,3) 
consider  the  designs 


D, 


C  0  0 
/  0  0 
U  2 


0  0 
0  0 
3  4 


111 

2  2  3 

3  4  4 


(  o  o 

D2  =  <  0  0 

U  2 


0  0  0  0 
0  0  0  0 
3  4  12 


0 

0 

3 


with  XQ  =2,  ^  =  2,  Aq  ^  =  4,  A^  =  0;  designs  and  D 2 

are  each  unions  of  generator  designs,  and  it  is  easy  to  verify  that  I>2 
is  inadmissible  with  respect  to  D^.  However,  D0  u  is  admissible 
with  respect  to  u  D3  where 


D 


3 


/‘  ‘  ‘ 
(223 

L  3  4  4 


and  hence  is  not  S-inadmissible  with  respect  to  D  .  A  sufficient 

condition  for  a  design  to  be  S-inadmissible  with  respect  to  a  design 

is  that  Aq  =  Aq  and  A^  ^  >  A^  ^ .  (See  Remark  5.1  of  [3J.) 

Thus  for  the  generator  designs 


8 


with 

that 


(1) 


AU)  =  i 
A1 


X(2)  -*  2 
0  ' 


^2)  =  o. 


D,  is  S-inadmissibie  with  respect  to  D  . 


respectively ,  we  see 


3.2  Exact  optimization  problem 

We  now  assume  that  for  given  (p,k)  there  are  n  generator  designs 

D.  (1  <  i  <  n)  no  two  of  which  are  equivalent,  and  no  one  of  which  is 

equivalent  to  the  union  of  two  or  more  generator  designs:  we  further 

assume  that  no  D.  (1  <  i  <  n)  is  S- inadmissible .  Let  A^\ 

be  the  design  parameters  associated  with  D^,  and  let  b^  be  the 

number  of  blocks  required  by  D.  (1  <  i  <  n).  Then  a  BTIB  design 
n  1 

D  =  u  f.D.  obtained  by  forming  unions  of  f,  >  0  replications  of 
i=l  n  ,  .  .  n  . . 

has  the  design  parameters  A^  =  £  f^O  anc*  A1  =  ^  ^i 'l* 

n  i  =  l  i=l 

and  requires  b  =  £  f.b.  blocks.  We  shall  consider  only  imple- 

i=l 

men table  D,  i.e.,  those  for  which  A^  >  0.  It  should  be  noted  that 

for  given  ^  =  (D^,...,D  )  the  design  D  is  completely  determined  by 

its  frequency  vector  £  =  (f  ,...,f  ). 

In  this  setup  the  integral  expression  (2.5)  for  the  confidence 

coefficient  can  be  regarded  as  a  function  of  £  for  given 
n 

p,  k,  b  =  £  f.b.  and  for  specified  £  =  dt/kb/o.  We  therefore 

i  =  l  11 

denote  (2.7)  by  g^.(^|^;  p,k,b;  £)  =  gr  (say)  where  the  subscript  E 
stands  for  exact.  Our  optimization  problem  can  be  stated  as: 


For  given  (p,k,b),  generator  designs  and  specified 


£  1  d*kb  o  >  0  choose  £  so  as  to 


maximize  g  (f|D;  p,k,b;  O 


subject  to  £  f.b^  =  b. 


i  =  l 


u  /  .  . 

.1  b'i  ’  »• 

i  =  l 


f .  >  0  (1  ^  i  <  n), 


(3.1) 


Sometimes  a  dual  of  (3.1)  might  be  of  interest.  This  happens  when 
b  is  at  the  disposal  of  the  experimenter,  and  he  wishes  to  guarantee  a 
specified  confidence  coefficient  using  the  smallest  possible  b.  In 
this  situation  our  optimization  problem  can  be  stated  as: 


For  given  (p,k),  generator  designs  and  specified 

d/o  >  0  and  a  (0  <  a  <  1)  choose  f  so  as  to 


n 

minimize  b  =  T  f.b. 

i  =  l  1  1 


subject  to  g  (£|j);  p,k,b;  £)  _>  1-a . 


.  L .  l  0 
1  =  1 


( 3 . : ) 


f.  '0  ( 1  <  i  <  n ) . 

^  -  .  r*  -Ml 


We  provide  numerical  solutions  to  (3.1)  for  selected  (p,k,b),  d/a 


3.3  Simplification  of  the  optimization  problem  for 


In  order  to  solve  (3.1)  it  is  first  necessary  to  list  all  generator 


D  for  given  (p,k)  and  then  to  enumerate  all  fre 


f  )  associated  with  admissible  designs 


D  =  u  f.D.  for  given  b  =  ) 

i=l  1  1  i=l 

such  as  o  =  4,  k  =  3  there  are 


Even  for  small  values  of 


manv 


and  the  number  n  of  such  designs  is  presently  unknown  to  us.  (It  is 


known  that  n 


different  methods  of  constructing  generator  designs  (these  being 


described  in  Section  3.2  of  [3])  but  it  is  known  that  such  methods  do  not 


it  turns  out 


each  such 


this  is  obviou, 


since  the  onlv 


two  generator  designs  possible  are 


(In  (3.3)  and  (3.4),  and  from  Section  3.3  on,  we  use  the  notation  D 


to  denote  a  generator  design  containing  zeros,  and  D  to  denote  a 


we  must 


consider  only  the  two  generator  designs 


0  0  0 


Do  ‘  \  1  1  2 

2  3  3 


(3.4) 


all  other  generator  designs  are  either  S-inadmissible  w.r.t.  D  for 

1  J  0 

b  =  3,  or  are  equivalent  to  u  f.D.  for  b  =  7  f  b  >3  For 

i=0  11  i=0  1  1 

example 


,(D 


0  0  0 
0  0  0 
12  3 


,(2) 


^3) 


are  generator  designs  which  are  S-inadmissible  w.r.t.  D  .  Also 


,(4) 


,(5) 


(3.6) 


are  generator  designs  which  are  equivalent  to  DQ  u  2D  for  b  =  5; 
therefore  any  BTIB  design  containing  D^4)  or  D^5)  (or  any  other 
equivalent  design  for  b  =  5)  can  be  replaced  by  an  equivalent  design 
employing  D,  and  Dj  in  place  of  D<4’  or  D<5>.  Thus  „e  can  limit 
consideration  to  D0  and  ^  of  (3.4)  to  construct  the  admissible 
designs  for  p  =  3,  k  =  3. 


Tor  given  (p,k)  suppose  that  we  have  only  two  generator  designs 


Dj ,  to  consider  as  in  (3.3)  and  (3.4)  for  p  >  2,  k  =  2  and 

p  =  3,  k  =  3,  respectively.  Let  X^\  t*le  Parameters 

associated  with  D. ,  and  let  b.  be  the  number  of  blocks  required 
t  l 

by  D.  (i  =  1,2);  note  that  =  0.  For  fixed  b  consider  all 

l  0  -  ^ 

possible  BTIB  designs  D  =  fQD0  u  f^D^  (f  ^  1)  with  £  f.b^  = 

i=0 

For  such  D  we  now  state  a  general  theorem  which  enables  us  to  rule 
out  as  inadmissible  certain  designs  with  sufficiently  "large"  f  ^ 

O 

("small"  f ^ ) .  We  first  note  that  for  such  D  we  can  express  n  and 
p  as 


2  .  .  ,(0)  (0)  (1)  , .(1), 

2  k  b{f0{bL\0  4b1X1  -bQX1  )  +  bX1  } 

n  U°  "  P  fK  i(0)  *  r k  J0)  K  .(i)^  A  K^(1)\ 

f0X0  {f0(blX0  +  pCblXl  ’Vi  ])  +  pbXl  } 


(3.7) 


P(f0} 


fo<bi‘[0>-bo4l)>  *  “l1’ 

wr  vr-vi1’’  *  -'i11 


(3.8) 


Theorem  3.1:  Let  (p,k,b)  and  (Dq,D^)  be  given.  For  these 
define 


c=  <p-DvS0>xi1> 


■  bi(xo°X0))(xo0)+xi0)>- 


(3.9) 


Consider  the  class  of  all  designs  D  =  fQDQ  u  f^D^  (f^  >  1,  >  0) 

r  U 

satisfying  l  f.b^  =  b.  Let  f  denote  the  upper  bound  on  fQ.  Then 
i=0 

we  have  either  of  the  following  two  cases: 
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Case  I :  If  C  >  0  then  there  exists  an  integer  f“  £  2  which  is  the 

2  2 

smallest  value  of  fQ  satisfying  n  (fQ)  >,  n  (fQ-a)  where  £  is 

the  smallest  positive  integer  such  that  divides  (abQ).  If 

f  <  f ^  then  all  designs  D  with  f,  >  f ^  are  inadmissible;  i.e., 

J  u  o  ***  0  - 

f“  -  1  is  the  largest  admissible  value  of  f  . 


Case  2 :  If  C  <  0  then  all  designs  D  are  admissible. 

The  above  result  follows  from  the  fact  that  n~  is  a  quasiconvex 
function  of  f  ^ ,  and  p  is  a  strictly  decreasing  function  of  f  . 

The  proof  of  the  theorem  is  given  in  Appendix  1.  If  Case  1  holds, 
in  order  that  f ^  <  f1"'  it  is  necessarv  that  b  be  sufficientlv 
large;  as  b  increases,  the  number  of  designs  eliminated  as  being 
inadmissible  increases. 


4.  EXACT  (DISCRETE)  OFTIMAL  DESIGNS 

4.1  General  results  for  p  i  2,  k  =  2 

For  p  ^  2,  k  =  2  any  BTIB  design  D  can  be  written  as  f^D  u  f  D 
where  f^  >  1,  f^>  0  and  are  given  by  (3.3).  For  we 

have  ^  =  1,  =  0,  bQ  =  p;  for  we  have  =  =  1 

b^  =  p(p-l)/2.  Thus  (see  (3.9))  C  =  p(p-l)/2  >  0,  and  Case  1  holds  for 
all  ?  >  2. 

General  expressions  for  n*’  and  p  are 


4b{2b  +  p ( p—  3 ) f _ } 
0 


0'  pf . { 2b  -  (p+l)f  ) 

u  u 


_  %  _  2(b-pfQ ) 
2b  +  p(p-3)f 


0 


(4.1a) 


• ■»  OV; 
ttW.r 


(  4  .  lb  ) 


is  the  smallest  fn  satisfying  (4.1c)  for 


'he  critical  number 


where  a 


is  odd:  thus  we  have 


where  in 


We  now  consider  in  detail  the  special  case 


to  illustrate  our  approach.  The  same  approach  is  used  for  the  othe 


D,  (of  (3.3))  for  arbitrary  b 


where  1  <  f .  <  b/ 


even 


according  as  b  is  even  or  odd,  all  inadmiss ible  f 


Thus  for  b  <  5  all  f,  are  admissible 


is  inadmissible;  for  b 
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and  (10,0)  are  inadmissible ,  etc.  In  Table  4.1  we  have  enumerated  all 

•  2 
designs  for  b  =  2(1)13,  and  have  given  the  associated  n  and  p; 

the  inadmissible  designs  are  noted. 

Conceptually  one  then  proceeds  as  follows:  We  are  given  (p,k,b), 
and  d/o  is  specified.  For  the  given  b  we  prepare  a  table  of  all  BTIB 
designs,  and  for  each  such  admissible  design  compute  the  associated 
confidence  coefficient  as  a  function  of  d/o.  This  has  been  done  for 
b  =  13  in  Table  4.2.  The  optimal  design  is  then  the  one  which  maxi¬ 
mizes  the  confidence  coefficient  for  the  specified  d/o.  We  note  in 
Table  4.2  that  each  admissible  design  is  optimal  for  one  of  the  tabulated 
values  of  d/o.  (We  do  not  claim  that  it  is  always  true  for  arbitrary 
(p,k,b)  that  each  admissible  design  is  optimal  for  some  d/o.)  Note 

that  the  design  f.  =  6,  f,  =  1  is  inadmissible.  Thus  we  have 
C  1 

obtained  a  numerical  solution  to  (3.1)  for  p  =  2,  k  =  2,  b  =  13. 

The  relevant  information  in  Table  4.2  concerning  the  optimal  designs  can  be 
summarized  as  in  Table  4.3;  in  this  table  the  optimal  f ^  are  denoted  by 
f.  (i  =  0,1)  and  the  associated  maximum  confidence  coefficient  by  g  . 
Finally,  tables  such  as  Table  4.2  can  be  prepared  for  arbitrary  b  ^  2 , 
and  each  of  these  can  be  summarized  as  in  Table  4.3.  Table  A3.1-I  in 
Appendix  3  combines  such  tables  for  b  >.  2  and  lists  for  each  b  =  2(1)19 
and  specified  d/o  =  0.2(0. 2)2.0  and  for  each  b  =  20(1)40  and 
d/o  =  0.1(0. 1)1.0  the  optimal  design  and  associated  g  . 

If  it  is  desired  to  solve  the  dual  problem  (3.2),  then  the  informa¬ 
tion  in  Table  A3.1-I  can  be  used  to  obtain  the  optimal  design  (minimum  b) 
to  achieve  a  specified  confidence  coefficient,  provided  that  b  <  40 
solves  the  problem.  If  b  >  40  is  required,  then  the  optimal  design 
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Table  4.2 

Confidence  Coefficient—  as  a  function  of  fA  and  d/o 
for  p  =  2 ,  k  =  2  when  b  =  13 


d 

a 

— 

6  = 

0 

— 

£  —  1 
o  ■  A 

ro  '  ‘ 

fo  =  3 

fli 

0  ’  ~ 

2 

n  =  27.13 

p  =  0.917 

n\  =  14. 3C 

o  =  0.913 

n"  =  10.20 

c  =  0 . 700 

n  =  9 . 36 

o  =  0.566 

rf  =  7.56 

p  =0.775 

H 

0.05 

0.255 

0.4542M 

0.4296 

0.4057 

0.3795 

0.3490 

0.3109  H 

0.10 

0.510 

0.4739M 

0.4571 

0.4385 

0.4164 

0.3834 

0.3509 

0.15 

0.765 

0.4936M 

0.4847 

0.4719 

0.4539 

0.4288 

0.3926 

0.20 

1.020 

0.5134M 

0.5124 

0.5053 

0.4918 

0.4700 

0.4354 

0.25 

1.275 

0.5332 

0. 5400M 

0.5383 

0.5298 

0.5114 

0.4737 

0.30 

1.530 

0.5529 

0.5675 

0.5719M 

0.5675 

0.5526 

0.5222 

0.35 

1.735 

0.5724 

0.5946 

0.6045M 

0.6045 

0.5931 

0.5653  '  [ 

0.40 

2.040 

0.5913 

0.6213 

0.6364 

0.6406M 

0.6326 

0.6074 

0.45 

2 . 295 

0.6109 

0 . 6474 

0 . 6674 

0.6754M 

0.6707 

0.6482 

0.50 

2.550 

0.6293 

0.6729 

0.6973 

0 . 7087M 

0.7070 

0.6871 

0.55 

2.304 

0.6484 

0.6974 

0.7259 

0.7403 

0.7413M 

0.7240 

0.60 

3.059 

0.6666 

0.7212 

0.7530 

0.7700 

0.7733M 

0.7584 

0.65 

3.314 

0.6844 

0.7440 

0.7787 

0.7977 

0.3028M 

0.7902 

0.70 

3.569 

0.7018 

0.7659 

0.3027 

0.8233 

0.3299M 

0.8193 

0.75 

3.824 

0.7188 

0.7866 

0.3251 

0.8466 

0.8544M 

0 . 8457 

11 

0.30 

4.079 

0.7353 

0.8062 

0 . 3458 

0.3679 

0.9764M 

0.S692 

0.85 

4.334 

0.7513 

0.8246 

0.3648 

0.8869 

0.8959M 

0.3901 

0.90 

4.539 

0.7667 

0.8420 

0.3821 

0.9040 

0.9130M 

0.9084 

0.35 

4 . 844 

0.7816 

0.3581 

0.8978 

0.9190 

0.9279M 

0.9242 

1.00 

5.099 

0.7959 

0.8731 

0.9119 

0.9322 

0.9407M 

0.9379 

—The  maximum  confidence  coefficient  in  each  row  is  marked  with  an  M;  the 

associated  f  and  f.  =(13-2f  )/2  denoted  bv  (f.,f, ),  is  the  ortimal  de- 
0  1  0  0  1 


a?o-n  fnr*  that  value  of  d/o  when  h  =  1.1. 
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Table  4.3 

Optimal  Design  and  Associated  Confidence  Coefficient 
as  a  Function  d/a 
for  p  =  2,  k  =  2  when  b  =  13 


r  o  o  " 

V'il  2. 

)•  *< 

13  -  2f  *  f 

0  1 

d/a 

fo 

0.05 

1 

1 

11 

0.4542 

0.10 

1 

1 

ii 

0.4739 

0.15 

1 

ii 

0.4936 

0.20 

1 

ii 

0.5134 

0.25 

2 

9 

0.5400 

0.30 

3 

7 

0.5719 

0.35 

3 

7 

0.6045 

0.40 

4 

5 

0.6406 

0.45 

4 

5 

0.6754 

0.50 

i 

4 

5 

0.7087 

0.55 

5 

3 

0.7413 

0.60 

5 

3 

0.7733 

0.65 

5 

3 

0.3028 

0.70 

5 

3 

0.8299 

0.75 

5 

3 

0.3544 

0.80 

5 

3 

0.8764 

0. 85 

5 

3 

0.8959 

0.90 

5 

3 

0.9130 

0.95 

5 

3 

0.9279 

1.00  j 

5 

3 

0.9407 
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for  1-a  =  0.80(0.05)0.95,0.99  and  d/o  =  0.2(0. 2)2.0  car.  be  found 
in  Table  A3.1-II.  (The  entries  in  this  table  were  obtained  by  a 
complete  computer  search  for  the  optimal  b.) 

It  might  be  noted  that  in  Table  A3.1-I  the  confidence  coefficient 
is  strictly  increasing  in  b  for  all  values  of  d/o.  That  this  is  not 
the  case  for  all  (p,k)  can  be  seen,  e.g.,  for  p  =  4,  k  =  2  in 
Table  A3. 3-1. 

4.1.2  Results  for  3<.p<6,k  =  2 

The  same  approach  as  was  adopted  to  obtain  Tables  A3.1-I  and 
A3.1-II  for  p  =  2,  k  =  2  was  also  used  to  obtain  corresponding  tables 
for  p  =  3(1)6,  k  =  2.  The  range  of  b-  and  d/o- values  covered 
in  each  of  these  tables  is  summarized  for  easy  reference  in  Table  4.4 
below. 


Table  4.4 

List  of  Tables  of  Exact  (Discrete)  Optimal  Designs  for  k  =  2 


p 

Optimal  Design 
and  g 

i - ] 

Optimal  Design 
to  achieve  1-a 

Range  of  b-,  d/o-  and  (l-a)-values 
in  each  table 

2 

s 

b  =  2(1)19  with  d/o  =  0.2(0. 2)2.0 
b  =  20(1)40  with  d/o  =  0.1(0. 1)1.0 

Table  A3.1-II 

1-a  =  0.90(0.05)0.95 ,0.99 
d/o  =  0.2(0. 2)2.0 

3 

Table  A3. 2-1 

b  =  3(3)46  with  d/o  =  0.2(0. 2)2.0 
b  =  51(3)114  with  d/o  =  0.1(0. 1)1.0 

1-a  =  0.30(0.05)0.95.0.99 
d/0  =  0.2(0. 2)2.0 

4 

Table  A3. 3-1 

D  =  **,3(2)36  with  d/o  =  0.2(0. 2)2.0 
b  =  36(2)76  with  d/o  =  0.1(0. 1)1.0 

Table  A3.3-II 

i-a  =  0. 30(0. 05)0. 95 ,0.99 
d/o  =  0.2(0. 2)2.0 

5 

Table  A3. 4-1 

b  =  5(5)85  with  d/o  =  0.2(0. 2)2.0 
b  =  90(5)200  with  d/o  =  0.1(0. 1)1.0 

Table  A3.4-II 

1-a  =  0.30(0.05)0.95,0.99 
d/a  =  0.2(0. 2)2.0 

6 

Table  A3. 5-1 

b  =  6^6)19(3)57  with  d/o  =  0.2<0. 2)2.0 
b  *  60(3)123  with  d/o  =  0.1(0. 1)1.0 

Table  A3.5-II 

1-a  =  0.60(0.05)0.95,0.99 
d/a  =  0.2(0. 2)2.0 

.  -  M _ i ,  ■  i  —  m 


jmtam 
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It  someriirses  occurs  that  g  decreases  as  b  increases  for  parti¬ 
cular  values  of  d/a  (see,  e.g..  Table  A3 . 3-1) ;  when  this  occurs  the 
experimenter  should,  of  course,  use  the  smallest  value  of  b  which 
guarantees  an  acceptable  value  of  1-a.  It  should  also  be  noted  that 
for  fixed  d/o  and  b  sufficiently  large,  the  confidence  coefficient 
gj.  is  strictly  increasing  in  b. 


4.2  Results  for  p  =  3,  k  =  3 

As  discussed  in  Section  3.3,  for  p  =  3,  k  =  3  we  can  limit  con¬ 
sideration  to  BTIB  designs  D  which  can  be  written  as  f  D  u  f,D 

CO  11 

where  f  >  1,  f  >  0  and  D  ,  D  are  given  by  (3.4)  for  arbitrary 

0  1  0  1 

b  =  3f  +  f,  (b  =  3,4,. . . );  here  1  <.  f  <  b/3,  0  <  f<  b-3.  For 

0  1  0 

Dg  we  have  =  2,  X^°^  =  l,  b^  =  3;  for  D,  we  have  A^J'^  =  0, 

=1,  b,  =  1.  Thus  (see  (3.9))  C  =  -3  <  0,  and  Case  2  holds. 
Therefore  all  designs  D  are  admissible. 

2 

For  p  =  3,  k  =  3  expressions  for  p  and  p  are 


2,  .  .  9b^ 

n  '  2f  ( 3b-4f  ) 

0  0 

b  -  2f„ 


P(f0)  = 


(4.3a) 


(4.3b) 


where 


b  =  3f0  +  fr 


(4.3c) 


Table  A3. 6-1  lists  for  each  b 


3(1)18  and  d/o  =  0.2(0. 2)2.0  and 


for  each  b  =  19(1)40  and  d/o  =  0.1(0. 1)1.0 


the  optimal  design  and 


associated  g^;  Table  A3.6-II  lists  the  optimal  design  for 
1-a  =  0.80(0.05)0.95,0.99  and  d/a  =  0.2(0. 2)2.0. 

5.  APPROXIMATE  (CONTINUOUS)  OPTIMAL  DESIGNS 


5.1  Preliminaries 

As  in  Section  3.3  we  continue  to  deal  with  situations  in  which  it 
suffices  to  consider  only  two  generator  designs.  In  Section  4  we  noted 
that  the  number  of  competing  admissible  designs  in  general  (but  not 
always)  increases  with  b  for  fixed  (p„k).  We  have  also  seen  that 
for  each  (p,k,b)  the  optimal  design  depends  on  d/a.  Thus  a  separate 
determination  of  the  optimal  design  must  be  made  tor  each  (p,k,b), 
d/a  combination.  This  represents  a  formidable  computing  and  tabulation 
task.  The  solutions  for  many  of  the  most  useful  combinations  are  given 
in  Appendix  3. 

In  order  to  extend  the  results  given  in  Appendix  3,  and  to  do  so 
in  a  compact  form  we  introduce  in  this  section  an  approximation  (for 
large  b)  to  the  discrete  optimisation  problem  considered  in  Section  4, 
i.e. ,  we  shall  replace  the  discrete  optimization  problems  (3.1)  and 
(3.2)  by  their  continuous  analogues.  The  continuous  versions  are 
analytically  more  tractable  and  computationally  easier  to  solve. 

Moreover,  since  their  solution  does  not  depend  individually  on  b 
and  d/o  but  only  on  these  quantities  through  $  =  d/kb/o ,  the  number 
of  solutions  that  must  be  tabulated  is  drastically  reduced.  The 
approximate  continuous  optimal  solutions  (discussed  in  Section  5.3  below) 
are  given  in  Tables  A4.1-I  through  A4.2-II;  in  Section  5.5  the  accuracy 
of  these  solutions  is  checked  against  the  corresponding  exact  discrete 
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w>  seek  tv'  obtain  >  as  thp  solution  in  >  ot  the  filiation  0, 

t;t  do  in*  so  kp  must  bp  insured  that  in  allowable  solution  in  >  exist  .  , 
that  it  is  unique,  and  that  it  is  in  tact  asaooiated  with  t ho  maximum  ,'t 
*,.•  Actually  it  turns  ovit  that  either  a  unique  solvit  ion  i;\  ^  >'t 

•'>  d  exists,  lies  in  the  interval  v  v' ,  1 ' ,  and  oonvsponds  tv<  the 


maximizing  value 


u'  that  the  maximum  value  ot 


for  >  t  [0,1] 


occurs  at  the  boundary  Y  =  0  for  f.  sufficiently  small  or  at  the 
boundary  >  =  1  for  (,  sufficiently  large;  the  solution  Y  z  1  occurs 
only  under'  Case  0  let.  Theoivm  3.1). 

We  show  in  Appendix  0  that 
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In  i 5 . 3 ) ,  fd  •  )  denotes  the  standard  normal  pdf  and  $^lx|t)  denotes  the 

c  it  at  the  equicoordinate  point  x  of  ait  r-variate  equicorrelated  standard 

normal  distribution  with  common  correlation  coefficient  t;  the  quantities 
■> 

n‘ .  o.  3n  /3y  and  3p/3>  are  given  as  functions  of  y  by  (A.D-1A.4), 
respectively,  in  Appendix  l.  Mote  that  the  sign  ot  3 g ^  3>  depends 
only  on  the  sign  of  h ( > | , D  ;  p  ,k ;  £). 

5.2.3  Study  of  g,  and  its  maximum 
_ A _ 

In  this  section  we  shall  study  the  behavior  of  g^  as  a  function  of 
Y  and  £  for  fixed  lp,k)  and  (fl^.D^).  It  is  straightforward  to 
check  that  in  the  limiting  case  >  =  0  we  have  n*  =  p  =  1  and 


hence 


rom 


All  of  our  calculations  and  studies  lead  us  to  conclude  that 


regarded  as  a  function  of  y  attains  a  unique  maximum  at  y,  the  value 


of  which  depends  on  £  and  (D  ,D  ).  For  all 


where  £  is  given  by  (5.6)  below,  we  have  that  g 


the  maximum  being 


equal  to  1/2.  This  result  parallels  those  obtained  in  [1]  and  [2] 


To  study  the  behavior  of  g  as  a  function  of  y  for  £  >  £ 


we  note  that  for  large  £  the  second  term  in  (5.3)  dominates  and  hence 


for  large  £  we  have  sgn(3g  /3y)  =  sgn(h)  =  -sgn( 3n  /3y).  In  Appendix  1 


we  have  shown  that  n 


is  a  quasiconvex  function  of  y  for  0  <  y  <.  1 


Depending  on  whether  n  achieves  a  minimum  in  the  interior  of  [0,1] 


or  at  the  boundary  y  -  1  we  have  the  following  two  cases  which  depend 


Case  1  ( C  >  0 ) :  In  this  case  g  has  a  unique  maximum  at  y  (0 


here  y  is  the  unique  solution  in  y  of  the  equation 


The  maximizing  solution  y  is  a  strictly  increasing  function  of  £  for 


In  the  limit  (£-*■«>)  the  maximizing  solution  approaches  y 


where 


which  is  the  largest  admissible  value  of  y.  This  limiting  value  y 
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also  be  easily  found  by  solving  the  quadratic  equation  3n  / =  0  in 
y,  i.e.,  g,  is  maximized  for  4  ■>  00  (which  causes  the  confidence 
coefficient  (1-a)  to  approach  unity)  by  minimizing  the  common  variance 
of  the  a.  -  a.  (1  «■  i  <  p). 

Case  2  (C  0):  In  this  case  there  exists  a  positive  constant 

5  =  ^1(DQ,Di;p,k),  (0<^0<C1<  ®>  such  that  for  all  4  <  Uq.^). 

g  has  a  unique  maximum  at  y  <  1;  here  y  is  the  unique  solution 
A 

in  y  to  (5.4).  The  maximizing  solution  y  is  a  strictly  increasing 

function  of  4  for  4  e  (4^.4^)  with  y  ■*  1  as  4  ■*  4^«  ('or  all 

(  >  4^  the  maximizing  solution  is  also  y  =  1  (which  implies  n£ 

replications  of  D^).  As  with  Case  1,  in  the  limit  (4  the 

maximizing  solution  is  the  value  of  y  which  minimizes  the  common 

variance  of  the  a„  -  a.  (1  <  i  <  p). 

0  1  “  ” 


5.2.4  Definition  of  4Q 

As  the  first  step  in  finding  a  closed  expression  for  4Q  we  consider 

lim  h ( y | D  , D  ;  p,k;  4).  From  (A. 3),  (A. 4)  and  (A. 5)  we  see  that 
Y-0  0  1  2 

lim  4^-  <  0  and  lim  4^—  =  -®  <  0.  Therefore 
Y-0  ^  ,*0  SV 


where  4Q  =  CQ( ;  p,k)  is  defined  by 


W<>  evaluate  this  limit  in  Appendix  A.'  and  how  it  to  bo 


t  plp-DK*  _  ,(0|  1 


b 


V  .pxJ 


1 

nr 


Note  that  E  depends  only  on  D  .  If  D  is  either  a  BIB  or  a  KB 
^  II 

design  between  the  p  tost  treatments  as  in  t3.3)  and  then 

i  i ) 

we  can  substitute  :  b,ktk-l)  ptp-l'  in  ts.6)  to  obtain 


1  „ 

-  —  p  at 

0  2  D  - 


toil  j )  ;  MfcAI- 

.  y  tk-lU 


t  r> 


Remark  3 .  I  :  It  is  known  ( see ,  e.g.,  Gupta  f5j)  that  *^0|l  at  - 
'fjtOii  3)  =  I  ,  *,(0|  1  3)  -  l  u  t  fl  :n  lares  inti  3),  and 

^  ,  k  0  |  1  J  >  =  1  '  8  +  ( 3  4n)arcsintl  3>;  $  tO|l.  3)  has  been  computed  to 

five  decimal  places  for  t  -  It  1)12  by  Gupta  [b,  Tab le  II,  p.  81']. 


B . 3 . >  Definition  of  £ 


For  Case  .  we  define  £ ,  as  the  smallest  value  of  £  tor  which 
>  equals  unity;  hence  for  £  -  ^  we  have  >  -  1 .  Thus  £  is  the 
solution  in  (.  of  the  equation 


h(y|  D0.P1 ;  p.k;  O 


'  0. 


I ti  general  f  depends  on  D  ,p  . 

l  0  1 

S.0.6  Uniqueness  of  maximum  of  g^  as  a  function  of  t 

As  mentioned  in  Section  b.3.3,  we  have  not  vet  proved  analvticallv 
the  existence  of  a  unique  maximum  for  g^  as  a  function  of  when 
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r  ^  <  r  v  ,0  (nor  vi id  we  prove  the  corresponding,  result  in  our  previous 

papers);  however,  all  of  our  numerical  calculations  and  certain  analytical 

considerations  point  to  this  conclusion. 

We  have  computed  g^  as  a  function  of  y  for  selected  values  of 

S,  and  given  the  results  in  Tables  5.1A  and  5. IB  to  illustrate  Cases  1 

and  2,  respectively.  Table  5.1A  is  for  p  =  k  =  2  (with  generator 

design;  (D,,B,)  of  (3.3)),  and  Table  5. IB  is  for  p  =  k  =  3  (with 
0  1  r 

generator  designs  (D^,D^)  of  (3.4));  these  computations  give  repre¬ 
sentative  pictures  of  the  behavior  of  g.  as  a  function  of  >  •  The 
behavior  of  g  in  Table  5.1A  is  analogous  to  that  of  g  in  Figure  1 
of  [l]  in  that  g^  has  a  unique  maximum  at  y  <■  y^  for  r  >  However, 

unlike  the  situation  in  Figure  1  where  lim  g  =  1/2P  we  now  have 

Y-l 

lim  g  depending  on  £  and  other'  parameters  of  the  design. 

_  i  A 

r*i 

5.3  Approximate  (continuous)  optimal  designs  for  given  (p,k),  (D^.D  ) 
and  specified  confidence  coefficient 

As  pointed  out  in  Section  3  (see  (3.2)),  sometimes  the  experimenter 
wishes  to  guarantee  a  specified  confidence  coefficient  1-a  for  given 
(p,k),  (D^,D^)  using  the  smallest  possible  b,  or  in  the  continuous 

case  using  the  smallest  possible  £,  sav  S.  This  dual  problem  can  be 
solved  in  the  continuous  case  in  the  following  wav: 

a)  If  Case  1  holds  or  if  Case  2  holds  but  £  <  £. ^  (see  comments 
below),  then  solve  simultaneously  the  two  equations 


h(Y|D0,D1;  p,k;  £)  =  0 


(5.8) 


CO 

r 

r  r-  L 

xv p  t  5 / n 

j  * 

CO 

_  fi-p 

dlK  x) 


1-a 


(5.8) 
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Values  of 


when  D 


0.7979 


as  a  Function  of 


when  D  =  f 

0  0 


4707 


5000 


4561 


4431 


4049 


6196 


6407 


6769 


6779 


6716 


6978 


7639 


8350 


0.0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

.5000 

.  4794 

.4630 

.4569 

.  4451 

.4321 

.4174 

.  4004 

.  3802 

.  3553 

.5000 

.5731 

.  5993 

.6161 

.6272 

.6334 

.6352 

.6321 

.6231 

.6063 
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for  y  and  £  ,  the  solutions  being  y  and  5 ;  here  h ,  n  and  p  are 

given  by  (5.3),  (A.l)  and  (A.2),  respectively. 

b)  If  Case  2  holds  and  (  >  then  solve  only  (5,9)  for  £  with 

Y  =  1,  the  solution  being  y  =  1,  £. 

To  see  whether  or  not  1  <  for  p  =  k  =  3,  (D^,D^)  given  by  (3.4), 

and  specified  1-a  can  be  found  by  consulting  Table  A4.2-I  in  the  Appendix. 

If  C  <  C.  for  the  smallest  tabulated  ^  1-a,  then  i  <  £  ;  if 

C  >  5^  for  the  largest  tabulated  gA  <_  1-a,  then  i  ^  5  .  For  the 

intermediate  case  in  which  the  S-values  corresponding  to  the  g  -values 

A 

bracketing  the  specified  1-a  are,  respectively,  <  5^  and  > 

then  it  is  necessary  to  use  trial  and  error  to  determine  if  £  <  or 

5.4  Use  of  tables  of  approximate  (continuous)  optimal  designs  for 
p  =  2(1)6,  k  =  2  and  p  =  3,  k  =  3 

For  given  (p,k)  and  (DQ,D^)  tables  of  approximate  continuous 
optimal  designs  can  be  computed  using  the  methods  described  in  Section  5.3. 
This  has  been  done  for  p  =  2(1)6,  k  =  2  for  (Dq,D^)  given  by  (3.3) 
and  for  p  =  k  =  3  for  (D^,D^)  given  by  (3.4).  The  results  are 
summarized  in  Tables  A4.1-I  and  A4.2-I,  respectively,  which  for  selected 
5  give  y  and  g  ,  and  in  Table  A4.1-II  and  A4.2-II  which  for 
1-a  =  0.75(0.05)0.95,  0.99  give  y  and  £. 

Tables  A4.1-I  and  A4.2-I  are  to  be  used  for  the  large  values  of  b  not 

.  2 
considered  m  Tables  A3.1-I  through  A3. 6-1.  Here  (p,k,b)  and  a  are 

given,  and  (p,k)  determines  (D^jD^).  The  experimenter  specifies  d 
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(his  "yardstick")  which  determines  £  =  dt/^/a.  Entering  the  appropriate 

table  with  (p,k,£)  the  experimenter  obtains  y  and  the  associated 

g  .  Then  using  the  relation  y  t  f  b  /b,  he  adopts  f  =  by/b  (to 
A  U  U  0  U 

the  nearest  integer)  as  his  (approximately  optimal)  design  which  will 

guarantee  a  confidence  coefficient  of  approximately  g  . 

A 

Tables  A4.1-II  andA4.2-II  are  to  be  used  for  the  large  values  of  b  not 

2 

considered  in  Tables  A3.1-II  through  A3.6-II.  Here  (p,k)  and  0“  are  given, 

and  (o,k)  determines  O  ,D The  experimenter  specifies  d  and  1-a. 

u  1 

Entering  the  table  with  (p,k,l-a)  the  experimenter  obtains  y  and  the 

.  , —  .  -.0 

associated  £.  Then  since  £  =  d/kb/o  he  computes  b  =  int[(£o/d)  /k]. 

Finally,  he  chooses  f„  so  that  b„f./b  -v  y  and  b„f„  +  b,f,  is  as 
close  as  possible  (dJ  b.  This  yields  an  approximately  optimal  design 
with  associated  confidence  coefficient  of  approximately  1-a. 

The  approximations  referred  to  in  the  paragraphs  above  arise  because 
we  use  a  discrete  design  which  is  as  "close  as  possible"  to  the  optimal 
continuous  design.  These  approximations  become  increasingly  more 
accurate  as  b  increases.  The  goodness  of  the  approximation  is  dis¬ 
cussed  in  the  next  section. 

5.5  Comparison  of  exact  and  approximate  optimal  designs 

To  illustrate  the  accuracy  of  the  approximation  provided  by  the 
continuous  optimal  designs  we  computed  the  exact  discrete  optimal  design 
and  the  approximate  continuous  optimal  design  for  p  =  2,  k  =  2,  d/o  =  0.2 
and  selected  values  of  b  (and  thus  O-  The  results  are  displayed  in 
Table  5.2.  In  this  table  we  have  denoted  b^^/b  by  y^,  where  f ^ 
corresponds  to  the  exact  optimal  design.  We  note  that  the  approximate 


able  5.: 


Comparison  of  Exact  and  Approximate  Optical  Designs 


(p  =  2,  k  =  2,  d/a  =  0.2) 


b 

1 

Exact  Optimal  Design 

Approximate 

Optimal  Design 

y'e 

Y 

eA 

10 

0.8944 

| 

0.2000 

0.5028 

■m 

0.5041 

12 

0 .  9~93 

fl 

0. 1667 

0.5104 

0.1737 

0.5104 

15 

1.0955 

fl 

0.2667 

0.5210 

| 

0.5210 

20 

1.2649 

fl 

0.4000 

0.5390 

0.3528 

0.5393 

25 

1.4142 

5 

0.4000 

0.5572 

0.4195 

0.5572 

30 

1.5492 

B 

0 . 4667 

0.5744 

0.4693 

0 . 5744 

35 

1.6733 

0.5143 

0.5083 

0.5908 

40 

1.7389 

B 

0.5500 

0.6063 

0.5398 

0 . 6064 

50 

2.0000 

15 

0.6000 

0.6352 

0. 5SS1 

0.6352 

75 

2.4495 

25 

0 . 6667 

0 . 6965 

0.6627 

0.6965 

100 


.  $234 


35 


0.7000 


0.7062 


3: 


continuous  optimal  designs  provide  excellent  approximations  to  the  exact 
optimal  designs  even  for  small  values  of  £  (associated  with  "low" 
joint  confidence  coefficients),  and  that  gr  is  only  slightly  smaller 
than  g^;  as  is  to  be  expected,  the  approximation  improves  with 
increasing 
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APPENDIX  1 


Proof  of  Theorem  3. 


sha 

the 

Sec 


Por  mathematical  convenience  and  without  loss  of  generality 
I  regard  y  -  f  ,b  b  as  a  continuous  variable  taking  values 
interva*  (3.1].  (We  use  the  same  continuous  approximation 


tion  5.)  for  p  '  k  =  2  and 


k  =  3  we  substituti 


b>  b,  in  1 1.’)  and  (3.8)  and  obtain 


n*(> )  1  —7 


k‘b  {y(b. \^01-b  \(0S  +  b^\(1)) 


*  b.\^'>]  ♦  ob  v!1'} 

11  C/-  1J  0  1 


ap.c 


0  (  >  >  = 


>ib  \j0).b  v5:>)  ♦  b  ' 

_ _ lx  J  1  U  1 


Hb 


b.\:° -b,\!i’.b.v,o’> .  b  a1'  ' 

a  1  Ox  X  0  01 


respectively.  It  follows  that 


3o 

3> 


b  b 

0  10  1 


{^b.Vp’-b  ♦  b  v^V  '  C 

11  01  10  01 


(1) 


for  b. .  \  s  0 ,  and  therefore  o  is  strictly  decreasing  in 


land  hence  in  f^).  Next  we  have 


3n~ 

3\ 


k‘b  iK>  3 


v^\‘{>CpCb1vi^-b0x^1>3  +  b^J01]  ♦  pb^:V 


where 
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,  v  .  2..  .(1)  ,  .(0)  ,  ,(0).r  ,,  .(0)  ,  .(1).  .  .(0), 

v(y)  -  Y  (bgi^  ”^1^1  -D1^0  ^1^1  "^0^1  +  blXQ 


T  V  JDr  I  K  K  *U-K 

2>Vl  [p(bixi  "Vi  }  + 


,  .(OK  ..  .Cl)  .2 
bixo  ]  -  p(boxi  )  ’ 


(A.  5) 


,  .  2  .  2 

Since  lim  n  =  ®,  it  follows  that  n  must  be  decreasing,  at  least 

Y-*0 

in  a  small  neighborhood  of  y  =  0+ ;  thus  it  suffices  to  show  that  n~ 

has  at  most  one  stationary  point  in  (0,1],  i.e.,  that  the  equation 

My)  =  0  has  at  most  one  root  in  (0,1]. 

(1)  2 

Since  the  constant  term  -p(b  X  )  in  (A. 5)  is  negative,  a 
necessary  (but  clearly  not  sufficient)  condition  for  both  roots  of  My) 
to  be  real,  positive  and  distinct  is  that  the  coefficients  of  y^  and  y 
in  My)  he  negative  and  positive,  respectively,  i.e.,  we  require  that 


P(b.x; 


(0) 


Vi  ) 


blX0 


<  0. 


Therefore  , 


3My ) 


3Y 


Y=1 


L(^°)  +  X[0))Cp(biXi0)-b0^1)) 


+  b 


*i0)] 


>  0. 


The  latter  inequality  shows  that  M y)  is  increasing  at  y  =  1.  This 

together  with  the  fact  that  My)  as  Y  -*■  ±°°  implies  that  at 

most  one  root  of  My)  must  be  in  (0,1]. 

The  proof  of  the  theorem  now  follows  easily  since  Case  1  or  Case  2 

holds  depending  on  whether 


sgn{b^C} 


sgn 


is  >0  or  <  0,  respectively,  where  C  is  given  by  (3.9).  If  Case  1 
2 

holds  then  n  first  decreases  and  then  increases  with  f^  (i.e.,  y) 

while  p  always  decreases  with  f  .  Hence  there  exists  a  critical 

L  o 

number  f^  which  is  the  smallest  value  of  fQ  at  which  n"  starts 

increasing  (for  fixed  b).  Thus  f^  is  the  smallest  value  of  f 
2  2 

satisfying  n  (fQ)  i  o~(fQ-a)  where  a  is  the  smallest  positive 
integer  such  that  ab^/b^  is  a  positive  integer.  If  f^  <  f ^ ,  then 
clearlv  designs  with  f„  >  f^  are  inadmissible.  If  Case  2  holds  then 

o 

since  both  n“  and  p  are  strictly  decreasing  in  fQ  (i.e.,  y), 

all  designs  D  =  u  f^D^  with  fQ  >  0  are  admissible. 


APFl'.NDI X 


Derivation  of  Results  in  beet  ion  '■ 


I.  Evaluation  and  simplification  of  .)> 

From  t  2 .  ’  1  and  the  definition  of  g.  g  i  von  in  Section  s  . .  1  we 

A 

obtain  bv  direct  calculation 


A* 

3y 


A  r  .p-1 

/  P* 


xnn'-K 
L  n* 1-p  J 


xn»  e  -t-c 

*1*1 -p  . 


ft  xl 


•1*  L-I'  vn  '  ro  t 


x*v 


•Ml 


x'uo+u(  n'  *  l-o  - 


n‘ ( l-o ) 


(  A .  6  1 


In  (,A.t>),  >ft  •  )  denotes  the  standard  normal  pdl  ,  fl  •  1  denotes  the 
standard  normal  cdf,  n'  -  An/Al  -  (1  2n)3n*  Ay  where  An"  Ay  is 
given  bv  ( A .  4 ) ,  and  p'  -  Ap  Ay  is  given  bv  (A.  3).  After  some  simpli¬ 
fication  l  A.  6)  can  be  written  as 


Ay 


f  2  «» 

/l_£l  f  x*r*L 

XTl*  p  K 

xn*p’+i. 

r 

-  n* 1-p  . 

V.  *p  -® 

_  n* 1-p  _ 

—  , _  “■ 

xn*  p  +  C 

_  n* 1-p  j 

-  n* 1 -p  _ 

f  V  X  1  d  X 


ft  x  Vix"'j 


ia.-m 


Making  the  change  of  variables 


V  *  (xnrp+C 1  n* l -p 


^  A . 7 )  can 


3y 


where 


and  ) 

bv  parts  in 


where 

From  (A. 12) 
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be  expressed  as 


—  E , 


2n  (1-p  )  ‘  K‘vp 


-[ 


S  n  p  ’  i 

tr  +  k  UVU-p^nr'}! 


R  >  p 


(A.  8) 


R  = 


1-p 


,  S  = 


n*l-p 


( A .  9  ) 


r  p- 1 

l  -  1  y*  (yH(y)$*(y)dy 


(a. io  y 


E2  3  /  (y)^(y )♦*( y)dy. 


(A. 11) 


denotf 


SS  ♦((y-S)  R).  We  now  evaluate  and  E , .  Integrating 


E1  with  U  =  L(y4*(y)  and  dV  =  y*(y)dv  we  obti 


sin 


r  -  1  _  S 

E1 - 3  Ei  +  —  E,  ♦  ( p-l)E 

R  R"  J 


( A.  12  ) 


r  p  -  2  2 

=  J  (y)<*  (y)$>'(v)dv 


(A. 13) 


we  have 


Ex  .  -JL  E„  *  E 

R‘>1  ‘  R'tl 


( A. 14) 


Substituting  (A. 14)  in  (A. 8)  we  obtain 
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3>  '  2n2(l-p2)3/2 


.  SnV 

^ptR2*! 


-  +  |  t2n'(l-p)-np'}  e  ♦  E 

)  R  J  ‘  r^(R‘  +  l) 


(A. 15) 


By  developments  similar  to  those  in  [1]  we  can  write 


2  t~2  \  J~2  J  P"1 

/R  +1  Vr+1' 


/(  2R2  +  1)(3R2  +  i”)  R  +1 


(A. 16) 


/(2R‘+l)2w 


2R2+1  /  P': 


'(2r2+i>(3r?7i)*3R%1 


(A. 17) 


Substituting  and  E^  from  (A. 16)  and  (A. 17)  in  (A. 15),  and  replacing 

R  and  S  by  their  definitions  in  (A. 9)  we  obtain 


__A 

9Y  '  , 


n‘  i' l-o 


-2Cn' yT^p*  ♦  <*>♦ 
n  p 


f  £  fPp  P  " 
-1  (_  n  v  1+p  1+p 


V’(p-l)  A  /  2  V  fi  /  1- 

'2tt(  ltpT  \n  V1+0/  p'2  ln 


-P  P 

( i+2p )  i+: 


h( Y I D  ,D  ;  P,k; 
2n 


where  h(y|Dg,D  ;  p,k;  £)  is  given  by  (5.3! 


II.  Evaluation  of  the  limit  in  expression  (5.5)  for  Eg- 
—  •» 

We  note  from  (A.l)  and  (A.  2)  that  lim  n*”  =  d0  and  im  p 

y-*0  y-*0 

Since  <J>  , ( 0 1 1/2 )  =  1/p  we  have  from  (5.5)  that 

P-1 


(  I  )  1 

Since  lim  (J/(y)  =  -pCb^A^  )  we  have 

Y+0 


1  a  l 

lim  -z— Y  2  =  k  —rrr  • 

/i-o2  i!L.  V2xi  p 


Substituting  (A. 22)  in  (A. 18)  we  obtain  the  desired  result  (5.6) 
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APPENDIX  3 

Tables  of  exact  (discrete)  optimal  designs 


Tables 

A3.1-I 

and 

A3.1-II 

(P 

= 

2, 

k 

S 

2) 

Tables 

A3. 2-1 

and 

A3.2-II 

(P 

s 

3, 

k 

= 

2) 

Tables 

A3. 3-1 

and 

A3.3-II 

(P 

= 

4, 

k 

= 

2) 

Tables 

A3. 4-1 

and 

A3.4-II 

(P 

r 

5, 

k 

s 

2) 

Tables 

A3. 5-1 

and 

A3.5-II 

(P 

= 

6, 

k 

= 

2) 

Tables 

A3. 6-1 

and 

A3.6-II 

(P 

= 

3, 

k 

= 

3) 

Table  A3.1-I 


Optimal  Design—  and  Associated  Confidence  Coefficient  gv 
as  a  function  of  d/o 
for  p  =  2,  k  =  2,  b  =  211)40 


■c  vQ. 


“r0  ‘1 


5 

0.6798 

3 

0.5319 

6 

0.6998 

3 

0.5356 

6 

0.6993 

.9683 


7  7  7  9 

0 . 8329  0.9256  0.9724  0.9915 


7  8  8  9 

0.9423  0.9327  0.9766  0.9933 


8  9  8 

0 , 9  o  "  9  o  .  o  o  o  T'  o  ,  j  o  o  9 


9  9  9 

0.9994  0.9997  1.0000 


1/ 

—  The  upper  entry  in  each  cell  is 


and  the  lower  entrv  is 


4 


'able  A3.1-I  (continued) 
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Table  A3.1-II 

Optimal  Design—  to  Achieve  a  Specified  Confidence  Coefficient 

as  a  Function  of  d/a 
for  p  =  2 ,  k  =  2 


b  =  «■  c 

D  ‘ro  ‘1 


T 


Confidence 

A  / 
u<  1 

0 

Coefficient 

(1-a) 

m 

B 

,  0.6 

DEI 

1.0 

J 

1.2 

1.4 

1.6  j 

1.8 

2.0 

616 

154 

69 

39 

13 

13 

■ 

3 

7 

0.99 

m 

41 

29 

16 

B 

8 

5 

B 

3 

3 

351 

S8 

39 

0 

15 

Q 

3 

6 

5 

4 

0.95 

144 

36 

16 

9 

6 

■ 

3 

0 

0 

241 

61 

27 

16 

wm 

n 

5 

5 

n 

3 

0.90 

96 

24 

11 

6 

fl 

D 

*■> 

B 

1 

179 

45 

20 

12 

3 

5 

5 

1 

3 

0.95 

69 

17 

8 

5 

3 

2 

B 

H 

i 

136 

34 

16 

9 

6 

5 

■ 

— 

| 

0.80 

51 

13 

6 

3 

2 

■ 

B 

B 

B 

X 

'able  A3. 2- 1 


Optimal  Design—  and  Associated  Confidence  Coefficient  gr 
as  a  Function  of  d/a 
for  p  =  3,  k  =  2,  b  =  3(3)114 


b  =  "3* 


3f 


1 


b 

d/a 

0.2 

0.4 

0.6 

1 _ _ _ 

0.8 

. 

1.0 

1.2 

1 . 4 

I  1.6 

| 

i  2.0 

1 

. 

' 

1 

m 

1 

1 

' 

1  1 

1 

1 

1  1 

■J 

0.17-21 

0.2285  | 

fmm 

0.3643 

0.4394 

1  0.5157 

0.5904 

0.6609 

0.7252 

0.7921 

X 

l 

i 

1 

|  i 

1 

1 

1 

1 

1 

0.3279 

0.4138 

0.5038 

!  0.5932 

0.6778 

(  0.7537 

0.8186 

0.8714 

0.9123 

0.9425 

9 

1 

0.3952 

"  1 

0.4770 

0.5700 

2 

0.6625 

2 

0.7650 

2 

0.8469 

2 

1  0.9066 

2 

0.9468 

2 

0.9717 

0.9859 

12 

1 

1 

3 

3 

3 

3 

3 

3 

0.4170 

0.5111 

0.6144 

0.7238 

0.3270 

0.9023 

0.9497 

0.9763 

0.9898 

0.9960 

15 

1 

2 

3 

3 

4 

4 

4 

4 

4 

0.4390 

0.5330 

0.6537 

i  0.7906 

0.3748 

0.9369 

0.9722 

0.9391 

0.9961 

0.9983 

19 

.  J 

1 

2 

3 

4 

4 

4 

4 

4 

4 

4 

0.4531 

0.5542 

0.6935 

1  0.8212 

0.9101 

0.9605 

0.9948 

0.9949 

0.9985 

0.9996 

o  i  i 

1 

3 

4 

5 

5 

5 

5 

5 

5 

5 

0 . 4646 

0.5741 

0.7273 

0.8537 

0.9347 

0.9751 

0.9919 

0.9977 

0.9995 

0.9999 

1 

3 

5 

5 

6 

6 

6 

6 

6 

6 

0.4738 

0.5954 

0.7563 

0.S804 

0.9521 

0.9841 

0.9956 

0.9990 

0.9998 

1.0000 

27 

1 

4 

6 

6 

7 

7 

7  I 

7 

7 

0.4813 

0.6154  I 

0.7923  | 

0.9030 

|  0.9647 

0.9397 

0.9976 

0.9995 

!  0.9999 

30 

\ 

X 

5 

6 

7 

7 

1 

/ 

7 

n 

7 

0.4376 

0.6340 

0.3067 

0.9210 

0.9744 

0.9934 

0.9986  1 

0.9998 

1.0000 

33  | 

1 

5 

7 

3 

8 

8 

8 

S 

|  0.4929  ! 

0.6519  1 

0.9290 

0.9354 

0.9813 

0.9958 

0.9993 

0.9999 

1 

36 

1  1 

6 

9 

9 

9 

a 

9 

a 

0.4976 

0.6696 

0.S469 

0.9470  | 

0.9363 

0.99T3 

0.9996  j 

1.0000 

39 

1 

*7 

9 

9 

10 

10 

10  ! 

0.5017 

|  0.6961 

0.8633 

0.9567 

0.9899 

0.99S2 

0.9998 

1 

uO 

10 

10 

10 

10 

10 

1 

0.5053 

1  0.7016 

0.3779  I 

0.9647 

0.9926 

0.9989 

0.9999 

4S 

2 

8 

10 

1  1 

X  X 

11 

11 

0.5092 

0.7163 

0.3913 

0.9711 

0.9946  j 

0.9993 

0.9OQQ  j 

48 

9 

11 

12 

12 

12 

12  ! 

0.5136 

0.7307 

0.9032 

0.9763 

0.9960 

0.9995 

1.0000 

1 

. 

1 - - 

1/ 

The  upper  entrv  in  each  cell  is  f  and  the  lower  entrv  is  g_. 

0  '  L 


Table  A3. 3-1 


i  continue..!'' 
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Table  A3.2-II 

Optimal  Design—  to  Achieve  a  Specified  Confidence  Coefficient 

as  a  Function  of  d/a 
for  p  =  3 ,  k  =  2 


/ 


f° 

0  0> 

1  2 ' 

\ 

°0 

2  3J 

•  °1 

=  {2 

3  3  . 

}• 

=  3fo 

+  3fl 

Confidence 

d/a 

Coefficient 

(1-a) 

0.2 

0.4 

0 . 6 

0.8 

1.0 

m 

a 

IB 

Q 

2.0 

0.99 

975 

246 

in 

63 

42 

30 

21 

n 

15 

12 

241 

61 

27 

16 

10 

7 

5 

1 

n 

4 

3 

0.95 

591 

150 

66 

39 

24 

mm 

m 

12 

9 

9 

141 

36 

16 

9 

6 

H 

H 

3 

2 

2 

0.90 

426 

108 

48 

27 

D 

12 

9 

5 

B 

6 

98 

25 

11 

6 

H 

3 

2 

2 

H 

1 

0.35 

330 

84 

39 

21 

15 

12 

9 

n 

6 

6 

73 

19 

9 

5 

3 

3 

2 

■ 

1 

1 

0.80 

264 

66 

30 

m 

12 

9 

6 

6 

6 

56 

14 

6 

■ 

3 

2 

1 

i 

1 

1 

i 

■1 


1/ 


—The  upper  entry  in  each  cell  is  b. 


and  the  lower  entry  is 


'O' 


4? 


Table  A3. 3-1 

Optimal  Designs—  and  Associated  Confidence  Coefficient  g„ 

as  a  Function  of  d/o 
for  p=4,k=2,b=  4,3(2)78 

P ,[° 0  3  °Y  .[  ■  ■  :  -  n. 

ll  2  3  4  j  1  1.2  3  4  3  4  4J 


! . 2312  0 . 3~5  2 


3 

0.1263 

0.2239 

3 

0.3493 

3 

0.4394 

3 

0.6265 

3  1 

1  0.7454 

3 

0. 33S1 

-*) 

0.2483 

0.3724 

0.5101 

0.6452 

0.7629 

HH 

0.9176 

0.3519 

jL 

C.4559 

1 

0.5623 

1 

0.6651 

1 

0.'562 

4 

0.S323 

4 

0.9079 

3  ! 

0 . 23^4 

m 

3 

0.5-90 

3 

0.7025 

3 

0.3245 

3 

0.9076 

3 

0. 9565 

0.3893  I  0.9303  |  0.958: 

0.971: 


1 

0.3913 

1 

0 . 4939 

3 

0.3039 

3 

0.4709 

0.6041  0.7350 


4 

0.7515 


3 


0.5103  0.6556  0.7929 


4  4 

0 . 6753  0.3260 


3 

0 . 6948 


2  5  3 

0.5446  0.7071  0.S57S 


4 

33  0.9657 


b  0 

0.7356  0.3334 


1  / 

—  The  upper  entrv  in  each  cell  is 


0 . 9112 


4 


0.3697  I  0.9406  0.9763 


3  3  3  3 

0.3335  0.9493  0.9300  0.9931 


0.9363  0.9961 


0 . 9994 


5 

0 . 0009 


0.9698  0.9902  0.9973  0.999*  0. 


b 

0.9926 


6 

0.9277 


5  5  5 

0.9434  0.9315  0.9950 


34 

1 

0.4353 

36 

3 

0.3319 

0.9459  0.9333 


0. 


and  the  lower  entry  is  g  . 
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Table  A3. 3-1  (continued) 
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3 
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Table  A3. 3- I I 


Optimal  Design—  to  Achieve  a  Specified  Confidence  Coefficient 
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The  upper  entry  in  each  cell  is  b. 


and  the  lower  entry  is 
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Table  3.4-1 
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Table  A3. 5-1 

Optimal  Design—^  and  Associated  Confidence  Coefficient 

as  a  Function  of  d/a 
for  p  =  6,  k  =  2,  b  =  6(6)18(3)123 
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Table  A3. 5- I I 

Optimal  Design  —  to  Achieve  a  Specified  Confidence  Coefficient 

as  a  Function  of  d/o 
for  p  =  6 ,  k  =  2 
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— ^The  upper  entry  in  each  cell  is  b,  and  the  lower  entry  is  f  . 
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Table  A3. 6-1 

1/ 
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—  The  upper  entry  in  each  cell  is  ^  and  the  lower  entry  is  g^,. 
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Optimal  Design—  to  Achieve  a  Specified  Confidence  Coefficient 

as  a  Function  of  d/a 
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—  The  upper  entry  in  each  cell  is  b,  and  the  lower  entry  is  f  . 


APPENDIX  4 

Tables  of  approximate  (continuous)  optimal  desi 

Tables  A4.1-I  and  A4.1-II  (p  =  2(1)6,  k  =  2 
Tables  A4.2-I  and  A4.2-II  (p  =  3,  k  =  3) 


Table  A4.1- 


"ontinuous  Optimal  Designs  as  a  Function  of  5 


for  k  =  2,  p  =  2(1)6 


(0  0 
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=  3 
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:  4 

r 

s 
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) .  7979 

50  =  ■ 

L.  6926 

*o  =  2-5214 

Y 

Y 

* 

SA 

Y 

®A 

1.00 

0.1395 

0.5121 

1.25 

0.3453 

0.5376 

1.50 

0.4522 

0.5681 

1.75 

0.5297 

0.6011 

0.0008 

0.5003 

2.00 

0.5381 

0.6352 

0.1618 

0.5168 

2.25 

0 . 6334 

0.6695 

0.2599 

0.5391 

2.50 

0.6694 

0.7032 

0.3375 

0.5650 

2.75 

0.6983 

0.7358 

0.4001 

0.5932 

0.0941 

0.5094 

3.00 

0.7219 

0.7668 

0.4514 

0.6226 

0.1782 

0.5273 

3.25 

0.7414 

0.7959 

0.4939 

0.6524 

0.2475 

0.5493 

3.50 

0.7576 

0.8228 

0.5296 

0.6822 

0.3053 

0.5738 

3.75 

0.7711 

0.8475 

0.5597 

0.7115 

0.3541 

0.5998 

4.00 

0.7825 

0.8699 

0.5854 

0.7399 

0.3956 

0.6267 

4.25 

0.7921 

0.3899 

0.6073 

0.7670 

0.4312 

0.6540 

4.50 

0.3002 

0.9076 

0.6262 

0.7927 

0.4619 

0.6812 

4.75 

0.8071 

0.9230 

0.6424 

0.8168 

0.4886 

0.7080 

5.00 

0.3129 

0.9365 

0.6565 

0.8391 

0.5118 

0.7341 

5.25 

0.8179 

0.9480 

0.6687 

0.8596 

0.5321 

0.7593 

5.50 

0.3222 

0.9577 

0.6794 

0.8783 

0.5499 

0.7833 

5.75 

0.3258 

0.9659 

0.6886 

0.8952 

0.5656 

0.8059 

6.00 

0.3239 

0.9728 

0 . 6967 

0.9103 

0.5794 

0.3272 

6.25 

0.3315 

0.9784 

0.7037 

0.9237 

0.5916 

0.8469 

6.50 

0.3337 

0.9830 

0.7099 

0.9355 

0.6024 

0.3651 

6.75 

0.8356 

0.9868 

0.7152 

0.9458 

0.6119 

0.8818 

7.00 

0.8372 

0.9898 

0.7199 

0.9548 

0.6204 

0.8970 

7.25 

0.3385 

0.9922 

0.7240 

0.9625 

0.6279 

0.9107 

7.50 

0.3397 

0.9940 

0.7276 

0.9690 

0.6345 

0.9229 

8.00 

0.8415 

0.9966 

0.7334 

0.9793 

0 . 6456 

0.9435 

8.50 

0.3427 

0.9981 

0.7378 

0.9865 

0.6543 

0.9594 

9.00 

0.8436 

0.9990 

0.7411 

0.9914 

0.6611 

0.3714 

10.00 

0.8446 

0.9998 

0.7454 

0.9967 

0.6705 

0.9866 

11.00 

0.8450 

0.9999 

0.7477 

0.9989 

0.6761 

0.9942 

12.00 

0.8452 

1.0000 

0.7489 

0.9996 

0.6794 

0.9976 

13.00 

0.8453 

1.0000 

0.7495 

0.9999 

0.6813 

0.9991 

00 

0.8453* 

1.0000 

0.7500* 

1.0000 

0.6833* 

1.0000 

*Largest  admissible  value  of  y  (i.e.,  yu). 


Table  A4.1-I  (continued) 
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3.50 
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0.1436 

0.5234 

4.00 

0.2040 

0.5432 

4.25 

0.2555 

0.5655 

0.0705 

0.5086 

4.50 

0.2997 

0 . 5894 

0.1319 

0.5239 

4.75 

0.3380 

0.6142 

0.1844 

0.5427 

5.00 

0.3714 

0.6396 

0.2297 

0.5637 

5.25 

0.4006 

0 . 6651 

0.2691 

0.5862 

5.50 

0.4262 

0.6905 

0.3037 

0 . 6096 

5.75 

0.4439 

0.7153 

0.3341 

0.6336 

6.00 

0.4690 

0.7395 

0.3610 

0.6578 

6.25 

0.4868 

0.7628 

0.3850 

0.6319 

6.50 

0.5026 

0.7851 

0.4063 

0.7056 

6.75 

0.5167 

0.3062 

0.4253 

0.7288 

7.00 

0.5293 

0.3261 

0.4424 

0.7513 

7.25 

0.5405 

0.8446 

0.4577 

0.7729 

7.50 

0.5506 

0.8619 

0.4714 

0.7935 

7.75 

0.5596 

0.8777 

0.4837 

0.8130 

8.00 

0.5676 

0.3923 

0.4949 

0.8314 

8.25 

0.5748 

0.9055 

0.5049 

0.3487 

8 . 50 

0.5813 

0.9175 

0.5139 

0.8647 

3.75 

l  0.5370 

0.9283 

0.5221 

0.3795 

9.00 

0.5922 

0.9379 

0.5294 

0.3931 

9.25 

0.5968 

0.9465 

0.5360 

0.9056 

9.50 

0.6010 

0.9541 

0.5420 

0.9169 

10.00 

0.6080 

0 . 9666 

0.5523 

0.9364 

10.50 

0.6135 

0.9761 

0 . 5606 

0.9520 

11.00 

0.6179 

0.9832 

0 . 5673 

0 . 9644 

12.00 

1  0.6242 

0.9920 

0.5772 

0.9312 

13.00 

0.6279 

0.9965 

0.5335 

0.9906 

14.00 

0.6302 

0.9985 

0.5374 

0.9955 

15.00 

0.6315 

0.9994 

0.5399 

0.9980 

16.00 

0.6322 

0.9998 

0.5913 

0.9991 

00 

0.6330* 

1.0000 

0.5932* 

1.0000 

’'‘Largest  admissible  value  of  y  (i.e.,  y^). 
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Table  A4.1-II 


Continuous  Optimal 


n  •  !/ 

Designs— 


to  Achieve  a  Specified  Confidence  Coefficient 


for  p  =  2(1)6,  k=2 


1 
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p-1 

P 


Confidence 
Coefficient 
(  1-0  ) 

P 

2 

3 

4 

5 

6 

7.0218 

8.8362 

10.3624 

11.7052 

12.9186 

0.99 

0.8373 

0.7401 

0.6729 

0.6226 

0.5831 

5.2989 

6.8621 

8.1885 

9.3623 

10.4276 

0.95 

0.8188 

0.7174 

0.6491 

0.5987 

0.5595 

4.3894 

5.8265 

7.0531 

3.1429 

9.1349 

0.90 

0.7968 

0.6912 

0.6220 

0.5718 

0.5331 

3.7764 

5.1219 

6.2909 

7.3259 

8.2700 

0.35 

0.7724 

0.6631 

0.5935 

0.5437 

0.5056 

3.2870 

4.5741 

5.6833 

6.6750 

7.5817 

0.80 

0.7440 

0.6312 

0.5616 

0.5126 

0.4756 

— ^The  upper  entry  in  each  cell  is  £,  and  the  lower  entry  is  y. 


Table  A4.2-I 


Continuous  Optimal  Designs  as  a  Function  of 
for  d  =  3 ,  k  =  3 
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1.0000 
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0.6441 
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5.25 
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0.9135 

8.00 
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0.6435 
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8.50 
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0.7973 
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5.75 
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0.9410 
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0.8401 

0.7118 

6.00 
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0.9518 

10.00 
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0.8839 

0.7444 

6.25 

1.0000 

0.9610 

11.00 

1.0000 

0.9205 

0.7754 

6.50 

1.0000 

0.9686 

0.9513 

0.8043 

6.75 

1.0000 

0.9750 

CO 

1.0000 

0.9802 
0. 9844 
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0.9906 
0.9928 
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0.9993 
0.999S 


Table  A4.2-II 


Continuous  Optimal  Designs-i  to  Achieve  a  Specified  Confidence  Coefficient 


for  p  =  3,  k  =  3 


r  o  o  r 
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